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Abstract. The Hadamard maximal determinant problem asks for the largest 
n X n determinant with entries ±1 . When n = 1 (mod 4) , the maximal excess 
construction of Farmakis & Kounias FK has been the most successful general 
method for constructing large (though seldom maximal) determinants. For cer- 
tain small n , however, still larger determinants have been known; several new 
records were recently reported in |OSDS| . Here, we define "3- normalized" n x n 
Hadamard matrices, and construct large-determinant matrices of order n + 1 
from them. Our constructions account for most of the previous "small n" records, 
and set new records when n — 37, 49, 65, 73, 77, 85, 93 , and 97, most of which 
are beyond the reach of the maximal excess technique. We conjecture that our 
n = 37 determinant, 72 x 9^^ x 2'^® , achieves the global maximum. 



Overview. How big can the determinant of an nxn real matrix be, given a bound 
on the size of its entries? This is the Hadamard maximal determinant problem, and 
an easy argument reduces it to that of maximizing the determinant on matrices 
with entries ±1 . From here on, we consider only such matrices^. Determinant- 
maximizing matrices of this type are important in the design of experiments; they 
are called a D-optimal designs. 

When n = 4:k + 1 , the determinant of a ±1 matrix can never exceed the bound 



due to Barba |Baj . In these orders, the maximal excess construction of Farmakis 
& Kounias |FKj has offered the most successful general method for construct- 
ing determinants achieving a large fraction of B(n) . For certain orders such as 
n = 13, 21, 25, 41, and 61, however, determinants beyond those attainable via 
maximal excess have been known for some time |R,a| l(]KM| IBh| IBHH| IBrj . In 
other dimensions, specifically n = 29, 33, 45, 53, 69, 77, 85 and 93 , a numerical 
ascent algorithm recently produced the largest known determinants |OSDSj . Here, 
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we present a general construction which accounts for many of the cases above. A 
sHght modification of our method accounts for the case n = 29 . Except for certain 
maximizers which have been obtained as incidence matrices for suitable block de- 
signs (n = 25 ,41 and 61), this leaves only the current n = 33 record unexplained. 
In addition to clarifying the structure of many current records, our construction im- 
proves the previous records in dimensions n = 77, 85 and 93, and sets new records 
in dimensions n = 49, 65, 73 , and 97. In a sequel to the present paper, we describe 
the modified construction which accounts for the case n = 29 , and show that for 
n = 37, it produces the determinant 72 x 9^"^ x 2'^^ = 0.945(37) . We conjecture 
that this value achieves the global maximum in its dimension. 

We base our construction on 3-normalized Hadamard matrices, which we define in 
^ below. Starting with a 3-normalized Hadamard matrix of order 4 A; , we make 
certain rank-1 modifications, and then construct a large-determinant matrix of order 
4A; -|- 1 by adjoining an additional row and column. 

To optimize our method in a given dimension, we seek a 3-normalized Hadamard 
"starting" matrix whose excess, defined as the sum of all matrix entries (see Defi- 
nition 1^21), is as large as possible. In ^ we derive an upper bound on this excess. 
The determinants attainable by the maximal excess construction are also bounded, 
given the bounds on the excess of an Hadamard matrix first noted by Brown & 
Spencer |BSj . and later improved by Kounias & Farmakis |K1<'| . Comparison of the 
analogous bounds suggests that our method has the potential to construct deter- 
minants larger than any obtainable using the maximal excess method in arbitrarily 
high dimensions. 

Furthermore, for n up to about 100, we find Hadamard matrices that attain or 
closely approach our bound, thus establishing the records mentioned above. On 
the other hand, we do not yet know of any infinite family of Hadamard matrices 
that attains our bound, while several such families are known for the maximal 
excess construction. We should also note that asymptotically, the largest fraction of 
Barba's bound that either method can attain is B{n)/V2 « 0.7lB{n) . Our bound 
does exceed that of Kounias & Farmakis by a lower order term as n — > oo , however, 
as shown in Proposition 14.41 

2. Preliminaries 

Maximal determinants. Let {±1}"^" denote the space of all nxn sign ma- 

2 

trices — matrices populated entirely by ±ls. This space has cardinality 2" < oo , 
so for each n S |Z| , there exists a sign matrix of maximal determinant. Define the 
maximal determinant function accordingly: 

(2) md:|Z|^|Z|, md(n) := max {det ^ : ^ € {il}"""" } . 
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The exact value of md(n) has been estabUshed for all n < 26 except n = 19, 22 
and 23. An elementary argument first given by Hadamard shows that we always 
have md(n) < nP'l'^ , with equality iff the maximizing matrix is Hadamard, i.e. a 
sign matrix with mutually orthogonal rows. A well-known conjecture states that 
Hadamard matrices of order n exist for all n = (mod 4) , so that the bound is 
achieved in these dimesions. This conjecture is known to hold for all n = (mod 4) 
with n < 428 , and many larger values as well |SYj . 

The present paper focuses entirely on the case n = 1 (mod 4) , for which md(n) 
remains largely unknown. For n > 29 with n = 1 (mod 4) , we only know md(n) 
when n takes the form + (a + 1)^ when a is an odd prime power !Brj. and 
when a = 4 • BHHj . For these values of n , as well as n = 1 , 5, 13, and 25, the 
maximizing matrix is essentially the incidence matrix of a suitable block design, and 
md(n) attains the upper bound B{n) for odd n given in above. Conversely, 
equality cannot hold here unless n = a^ + (a + l)^ for some integer a , which forces 
n = 1 (mod 4) . Conjecturally, Barba's bound is attained whenever n has this 
form. 



Tensor notation. Given vectors x, y € R'^ with coordinates Xj, yj , we will iden- 
tify the tensor product x y with the rank-1 linear transformation whose action 
on R*^ and matrix entries [x (X" y]ij are characterized respectively by 

(3) (x®y)(v) := (y • v)x , [x (g) y]ij = . 
For any A G Hom(R") , we then have 

(4) yl o (x (g) y) = (Ax) (g) y , and (x (g) y) o ^ = x (*Ay) , 

where *A denotes the transpose of A . In particular, when ^ = c (g d for some 
c, d G R" , we have 

(5) (x(gy) o (c d) = (y • c)x d . 
Note also that 

(6) *(x y) = y x . 

Using these facts, we get well-known formulae for the inverse and determinant of 
the "rank-1 update" of a non-singular matrix: 

Lemma 2.1 (Sherman-Morrison formulae). If x,y £ R"" and Anxn is invertible, 
then 

det(A + x(g)y) = (det^) (1 + A"^x • y) 
N-i .-1 A-ix®(U"i)y 
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Proof. First consider the special case A = In , and choose an orthonormal basis 
{sj} for R" with ei = x/|x| . Then 

(x (g) y){ej) = (y • ej)x = (y • e,-)|x|ei , 

It follows that, relative to this basis, the matrix representing Z„ + x y has zeros 
everywhere except along the first row, and the diagonal. In particular, it is upper 
triangular, with all diagonal entries equal to 1 except the first, which is 1 + x • y . 
The determinant of I„ + x (gi y then reduces to this single entry, and we then have, 
for any invertible n x n matrix A , 

det(^ + X y) = (det A) det + ^"^ (x ® y )) 
= (det A) det (/„ + (A~^x) ® y) 
= idetA){l + A-^x-y) . 

This gives the stated formula for determinants. 

One can verify the formula for the inverse directly: Multiply it on the right by (ylx 
y) , and simplify, using the rules in ^ The result, a straightforward calculation, is 

In. □ 

Definition 2.2. To streamline the display of vectors and matrices below, we in- 
troduce the following convention: Given any integer n, and any A; G |Z| , n^ := 
re(l, !,...,!) G R'^ . Abusing notation slightly, we also use n^ to indicate a sequence 
of k copies of the number n . For instance, 

(3ii, -Ig) = (3, 3, 3, -1, -1, -1) eR20 . 

^ V ' ^ V ' 

11 9 

Excess. If A £ {±1}"^"' , then the sum of the entries of A measures the extent to 
which the +ls "exceed" the number of —Is in A . 

Definition 2.3 (Excess and row sums). The excess of any vector or matrix X is 
the sum of its entries, and we denote it by ex (X) . Note that for a vector x G R'^ 
or an n x m matrix A respectively, we have 

(7) ex(x):=x-lfc, ex (A) := In ■ Aim ■ 

We call the excess of the ith row of a matrix its ith row sum. 

Observation 2.4. If :x. £ R" and y G R™ , then 

ex (x (g) y) = ex (x) ex (y) . 



This follows immediately from (jSj). 
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The injection E : R'^^" j^(n+i)x(n+i) _ Let R''^" denote the space of aU real 
nxn matrices. Several of the facts above come together in relation to the map 
E : R"^" j^(n+i)x(n+i) ^ -^Ljgj^ adds an additional first row and column to A by 
prepending —1 to each row of A, and then prepends an entire row of +ls to that. 
In block form, we may express this as follows: 



(8) EiAr 



1 1. 

-Ir, A 



Lemma 2.5. For any square matrix A , we have 

det(S(^)) = (detA) (l + ex(A-^)) . 

Proof. Adding the first row of E{A) to each successive row produces the matrix 

1 In 

without changing the determinant. The determinant of the latter matrix is simply 
the determinant of its lower-right block, however. We thus have 

det {E{A)) = det(A + 1„ 1„) . 

The desired result now follows immediately from Lemma l2.1l and (|7|). □ 



Hadamard matrices. As mentioned above, a Hadamard matrix is an element 
H G {±1}'^^" whose rows are mutually orthogonal. We denote the set of n x n 
Hadamard matrices by Had (n) . 

Row and column permutations, and row and column negations, clearly leave Had (n) 
invariant. Hence 

Definition 2.6. Hadamard matrices Hi, H2 G Had (n) are equivalent if Hi can 
be transformed into H2 by permuting and/or negating sets of rows and/or columns. 

We also note here that Had (n) C {±1}"^" is invariant under transposition. Indeed, 
the defining condition for H G Had (n) implies that H^H = nin , and hence that 
i^-i = tH/n . 

Proposition 2.7. If n > 2 and H G Had (n) , then n = 4k for some k G , 
and using only column swaps and column negations, we can put the first 3 rows of 
H into the following standard form 



(9) 
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Here '+' and ' respectively abbreviate and —Ik- Moreover, if we denote 

the ith row sum of H by ri , we have 

Ak 

(10) Y.'l = lQk^. 

1=1 



Proof. To begin, note that by negating column j ol H = [hij] whenever hij = — 1 , 
we get an equivalent matrix whose first row contains only +l's. If we then sort the 
columns of the resulting matrix in decreasing lexicographic order, left to right, the 
first three rows of the resulting ff-equivalent matrix will necessarily take the form 

+ lfci +lfc2 +1^3 +lfc4 \ 
+ — lfe2 +lfc3 — 1A:4 / 

where ki + k2 + k^ + ki = n , and the mutual orthogonality of these rows gives three 
additional equations: 



The unique solution of the resulting 4x4 system is A:i = A:2 = A:3 = A:4 = n/4 , and 
since each ki is an integer, we have n/4 G Z"*" . This proves the first two claims. 

To get the fact about row sums, note that the the vector (ri, r2, . . . , r^k) of row 
sums coincides with Hl4k ■ It follows that 

4fc 

^rf = HUk ■ HUk = Uk ■ {'HH)Uu = Uk ■ (4fcl4fc) = {^kf . 
1=1 



This completes the proof. 



ki + k2 — k^ — ki = 
ki — k2 + k^ — ki = 
ki — k2 — k^ + ki = 



□ 



Though the first three rows of any Hadamard matrix can always be put in the form 
©, a different normalization will better suit our purposes. 

Definition 2.8 (3- normalization). When the first three rows of a Hadamard matrix 
H have the form 




(11) 

and all row sums of H are non-negative, we say that H is 3-normalized. 
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Proposition 2.9. Every Hadamard matrix H £ Had (n) , n > 2 , can he 3- 

normalized using column swaps, column-negation, and row-negation only. Moroever, 
every 3-normalized Hadamard matrix has the form 

/ + - - + \ 

+ - + - 

+ + - - 

^ — a4 b4 C4 d4 

\ t)„ Cn dn J 

Here and '—' respectively abbreviate ±1^ , and a^, bj, Cj, dj E {±1}"/^, and 
satisfy, for each i = A, . . . ,n , 

(12) ex (aj) = ex (bj) = ex (cj) = ex (dj) . 

For the row sums ri , we have ri = r2 = r^ = , and, for i > 3 , ri = n (mod 8) . 
If n = A (mod 8) , r^ < n — 8 for all rows, unless n = A. 

Finally, we have ex (A^) = n (mod 8) , and when n = (mod 8) , ex (iV) = n 
(mod 16) . 

Proof. By Proposition 12.71 we can write n = Ak and put the first three rows of 
H into the form Q using only column swaps and negations. By further negating 
columns k -\- 1 through 3k , we obtain an i?-equivalent matrix whose first three 
rows have the form (jllj) . Negating all rows having negative row sums, we now get 
a 3-normalized Hadamard matrix equivalent to H. This proves the Lemma's first 
statement. 

To get (fT2]). define Oj = ex(aj) , 6j = ex(bj) , etc., and recall that for any fixed 
z = 4, 5, . . . , n , row i is orthogonal to each of the first three rows. Consequently, 

tti - hi - Ci -\- di = 
Oj - 6j + Cj - dj = 
o-i + hi — Ci — di = . 

The general solution here is ai = hi = Ci = di , which proves ()12() . It also shows 
that rj = (mod 4) for all i . In fact, the first three row sums clearly equal 0, and 
to get the more precise claim about rows 4 through n , simply note that since aj 
has k entries, each congruent to 1 mod 2, ex (aj) has the same parity as k does. 
It follows that Tj = 4ex (aj) = or 4 (mod 8) , depending on whether k is even or 
odd respectively. 

For odd /c > 1 , no row can consist entirely of Is (giving row sum n), because 
the mutual orthogonality of rows would then force all remaining row sums to equal 
0^4 (mod 8) . Therefore, no row sum exceeds n — 8 in this case, as claimed. 
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Since we have an odd number (n — 3) of non-zero row sums, all congruent to n 
(mod 8) , it follows that ex (N) = n (mod 8) . 

Finally, since ri = n (mod 8) , we may write = 8/i and n = 8m when n = 
(mod 8) . By definition of the the excess, and above, we then have 

1=1 i=l 

Since x'^ = x (mod 2) for any integer x , the two sums above, and m , all have the 
same parity. We conclude that ex (A^) = 8m (mod 16) , as claimed. 

□ 

Remark 2.10. When n > 20 and n = 4 (mod 8) , the proof of Lemma 2.2 in 
|Haj generalizes to show that the largest possible row sum allowed by Proposition 
12.91 above, namely n — 8 , can be achieved by at most one row. Modulo column 
negations and permutations, such a row forms, along with the three initial rows of 
the 3-normalized matrix to which it belongs, a Hall set, as defined by Kimura |Kij . 

3. A CONSTRUCTION 

Constructing r(A^) . We now describe the basic construction we use to produce 
large-determinant sign matrices of order n + 1 . 

Step 1: Select a 3-normalized Hadamard matrix G Had (n) , n = 4A; . 

Step 2: Alter by negating its first k columns to produce a new Hadamard 
matrix A^' . By Lemma 12.91 the first three rows of A^' will take the form 




Step 3: Make the rank-1 modification 

N' ^ N' + 2e3<S)ek=: N" , 

where, for any m S |Z| we define em ■= (Imj 0„_m) G • This modifies 
N' by changing the —Is in its upper left 3 x A; block to -|-ls. Note also that 
A^" is no longer Hadamard. 

Step 4: Construct the final sign matrix r(A^) of size n -|- 1 by applying the 
injection E defined in ^ 

r(A^) = E{N") . 
The following result suggests the efficacy of this construction. 
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Theorem 3.1. Let N be a 3-normalized Hadamard matrix of order n = 4k . Then 
the sign matrix T{N) of order n + 1 described above has determinant 

det T{N) = (det N) (2 + ^) = in"/^ (2 + ^) . 

Proof. Using the notation of our construction above, we can express detr(A^) in 
terms of N" by applying Lemma 12.51 to get: 

(14) detr(iV) = det TV" (l + ex {N"-^)) . 
We proceed to expand the right-hand side in terms of N itself. 

First, given the definition of N" in Step 3 of our construction, Lemma l2. II enables 
us to write 

(15) det A^" = (det N') (l + 2 N'~^e3 ■ e^) . 

Since N' is Hadamard, we have N' = W/n, and hence N' 63 = Wea/n. But 
We3 simply sums the the first three columns of W , ie. the first three rows of N' 
itself. Step 2 of our construction displays these rows explicitly, and we easily deduce 

(16) N'-\s = -{-3k, -l3fc) . 

n 

It follows that 

(17) l + 2iV'-^e3-e, = l + ^f )=l-^-3k = -l, 

n \ -Isfc y V ^3fe / n 2 

and consequently, 

(18) detiV" = (-l)'=+i^ . 

The factor (— l)'^ appears because we have replaced N' by , a matter of negating 
k columns. 

We next calculate ex {^N"~^) . 

To begin, invert N" = N' + 2e2, ® using Lemma l2. 11 and compute 

ex OV"-i = exfiV'-^) -2— ^ . L 

(19) = ex [n'-^) + 4ex (iV'"'e3) • ex {{'N'rht) • 

Here we used (|17|) to evaluate the denominator, and Observation 12.41 to factor the 
excess of the tensor product. 
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Now unpack the excess terms on the right above. For i = 1, 2, . . . , n , let rj and 
r- denote the ith row sums of and A^' respectively. Since A^' is Hadamard, we 
may then write 



1 1 



n ^ ' n 

i=l 



We obtained A^' from a 3-normalized Hadamard matrix A^ by negating columns 1 
through k. So by the definition of 3-normalized form, we have r[ = —2k = —n/2 
when z = 1,2,3. When i > 4 Proposition 12.91 similarlv implies r- = rj/2. Finally, 
since ri = r2 = rs = , we have Y17=4 '''i ~ ^11=1 ^« — (-^) • Thus 

_ 1 , V" _ 1 ( -^n + e^{N)\ _ ex(A^) 3 

2 ■ 



Given these facts, we now easily analyze ex ( A^' ^es) using ()16|1: 



(21) ,.(^,-,3).i(,_3)^.,-l)^) 



n , ^ , 3n \ 6 3 
~4 ~ ~2 



Lastly, we compute ex ((W') -"^efc) , recalling again that we produce A^' from A^ G 
Had (n) by negating columns 1 through k , so that 

i^N'y^ek = -N'ek = - ■ (sum of the first k columns of A^' ) 

n n 

= • (sum of the first k columns of A^ ) . 

n 

Since A^ is 3-normalized, rows 1, 2, and 3 have only -|-ls in these columns. The 
first three rows therefore contribute k = n/A each to the sum of columns above. 
Proposition 12.91 further guarantees that rows 4 through n each contribute exactly 
l/4th of their full row sum. Denoting the latter row sums by rj as above, we may 
therefore write 

{^N'Y^ek = -— (na, r^, rg, . . . , r„) , 
4n 

whence 

(22) ex {m-'eu) = -^^ (3 n + ex (AT)) = - \ , 

since Y11=a '''i ~ i^) ' noted in deriving (pn|) above. 

Finally, we assemble the various facts above to calculate det r(A^) . First, combine 
equations (UH), ((THl) and (UH), to get 

detr(Af) = (-l)^'+i(detA^)-^ - (l + ex^Af'"^) +4 ex (^N'^^s) • ex ((^A^O^^e^) 

Then replace the excess terms here with the results of equations 1)20(1 , l|21() , and (|22() . 
Routine simplification then verifies the theorem. □ 
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4. Upper bounds 

Compare the formula in Theorem 13.11 above with the analogous formula we get by 
applying Lemma 12.51 to an arbitrary n x n Hadamard matrix H : 

(23) det {E{H)) = (det H) (^1 + ^ ex (H)^ . 

The latter formula clearly shows that we maximize det E{H) over all H £ Had (n) 
by maximizing ex (H) . Indeed, this constitutes the maximal excess technique of 
Farmakis & Kounias IFKj we mention in our introduction. 

By comparison, Theorem 13.11 above shows that we can improve on that technique 
in any dimension n where there exists a 3-normalized n x n Hadamard matrix N 
such that 

(24) ex {N) + n> ex (H) for all H £ Had (n) . 

Of course, the 3-normalized matrices form a relatively small subset in any equiva- 
lence class of Hadamard matrices, and one expects the excess function to reach a 
smaller max on this subset than on the entire class. So it is not obvious that we can 
attain the desirable situation expressed by 

We have found, however, that 3-normalized Hadamard matrices satisfying (|24() do 
exist, at least in dimensions n < 100 . The excesses of arbitrary Hadamard matrices 
and 3-normalized Hadamard matrices satisfy certain natural upper bounds. We 
discuss both bounds below, and compare them. We find that for large n , the 
bounds themselves do satisfy (|24j) . offering the hope that one could find an actual 
3-normalized matrix that does so. 

In this section, we briefly review known bounds on the excess of a general Hadamard 
matrix H £ Had (n) , then derive a new upper bound for the excess of a 3-normalized 
Hadamard matrix. Our bound clarifies what we can (and cannot) hope to obtain 
from the construction in ^ and indicates how to optimize the search for suitable 
3-normalized Hadamard matrices, which we take up in ^ 

Proposition 14.41 at the end of this section gives our comparison of the two bounds. 

The n'^/^ bound on maximal excess. The excess of a general Hadamard matrix 
H £ Had (n) can never exceed n^/^ . Indeed, the Cauchy-Schwartz inequality gives 

ex (H) = In ■ Hln < |ln| • \Hln\ = \fn ■ \fn ■ a/u , 

since H/^/n is orthogonal. This fact was apparently first observed by Brown & 
Spencer ^BS^, and independently, by Best |Bej . The bound is sharp in the sense 
that regular Hadamard matrices actually attain this excess value, and such matrices 
occur in infinitely many orders (for instance, order whenever Had (n) is non- 
empty [HH])- 
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In |Ki^'j ■ Kounias & Farmakis derive an upper bound smaller than r?!'^ when n is 
not a perfect square, and tabulate their bound versus the largest known excess for 
all orders up to n = 100 |FKj . For convenience, we reproduce their list of bounds 
in Table 1, 33 

A 3-normalized excess bound. If is a 3-normalized Hadamard matrix of 
order n > 2 , Prop. 12.91 shows that the first three rows of N have row sum zero, 
while the remaining row sums belong to the set 

{{0,8, 16, ... ,n} , if n = (mod 8) , 

{4} , if n = 4, 

{4,12,20, ... ,n- 8} , if n = 4 (mod8),n>4. 



Observation 4.1. Consider the last n — 3 rows of a 3-normalized N E Had (n) . 
For each r G RS (n) , let n^. count how many of these rows have row sum r . Then 



E 

RS(n) 



n 



3 , rn^ 
RS(n) 



ex (A^) , and ^ 

RS(n) 



n 



Proof. The first identity is trivial. The second totals the row sums — the excesses — 
of the last n — 3 rows. Since the first three rows of N contribute zero excess, this 
gives ex (A^) . Equation (fT(i|) from Prop. [23 evaluates the third sum; we record it 
in slightly different notation here for convenience only. □ 

Notation. In the theorem below, the notations 

Ix] and [x] + 

respectively denote the ceiling and positive ceiling of x , ie. least integer, and the 
least positive integer, respectively, which exceeds or equals x . 

Theorem 4.2. Suppose N £ Had (n) is 3-normalized, with n > 4 . Then ex (N) < 
u* , where 



(25) 
and 



Pn ■-- 



8/ 



Pn{n-3) ^ (n-4)(n-12) 

2 2pn 

if n = (mod 
, if n = 4 (mod 



n 



n 



n 



Equality holds in \25\) iff the last n — 3 row sums of N take the values Pn^'^ only. 



Proof. Our main argument degenerates for the cases n = 4, 8 and 12; we conse- 
quently argue the latter separately after handling the generic case n > 12 . 
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The case n > 12. As in Observation 14.11 we let rir count how many, among the 
last n — 3 rows of , have row sum r for each r G RS (n) . We then adapt 
a technique appearing in Kounias and Farmakis KP^ by considering the following 
function F of p G RS (n) : 



(26) 



RS(n) 



r) > 



The inequality here is crucial. It holds because the are all non-negative, while 
{p — r) and (/5 + 8 — r) necessarily give consecutive integer multiples of 8, and hence 
a non-negative product. 



On the other hand, F{p) expands as 

F{p)= p{p + S)Y,nr-{2p + S)Y,rnr + Y.' 
and we can evaluate all three sums explicitly using Observation 14.11 After careful 
simplification, we may then rewrite the estimate H26|) above as 



(27) 



exiN)<Gnip) :-- 



{p + A){n-3) (n-4)(n-12) 



+ 



2 ' 2(p + 4) 

Note too, for later use, that equality obtains here iff F{p) = . 

To extract the strongest possible result from (|27|) . we now minimize Gn{p) over 
RS (n) . Here, the assumption n > 12 comes into play; one easily sees that in 
this case Gn is concave up, hence attains its minimum on RS (ti) at the smallest 
p £ RS (n) for which the forward difference G„(/9 + 8) — Gn{p) is non- negative. 
Accordingly, we compute 



G„(p + 8)-G„(p) = 4 



(p + 8)2(n-3) -n2 



(p + 4)(p + 12) 

and observe that the denominator here is positive for all p G RS (n) . This makes 
non-negativity of the forward difference on RS (n) equivalent to 

n 



P> 



\/n — 3 

Since RS (n) depends on the residue class of n mod 8, so does the minimal solution 
Pmin G RS (n) for this inequality, but a little thought shows that we can express it 
in the following way: 

y I- 

n 



Pn 



8 



8 



8/ 



n 



n 



if n = (mod 



8/ 



n 



if n 



(mod 



By evaluating with p = pmin and then substituting p„ 
obtain (|25() . the main result of our theorem. 



-|- 4 , we now 
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It remains to analyze the condition for equality there. According to the remark 
immediately following ^27\i . this condition corresponds to -F(pjnin) = 0, and by 
(|2())) . the latter will happen if and only if we have 



The factor (/Omin — r){pmm + 8 — r) here vanishes when and only when r = pmin 
or r = pmin + 8 , or equivalently, when r = p„ it 4 . Equality in consequently 
forces rir = for all other r G RS (n) , exactly as claimed. 

The cases n = 4, 8 and 12. Since RS (4) = RS (12) = {4}, the excess of a 3- 
normalized Hadamard matrix N will equal (4 — 3) x 4 = 4 if N £ Had (4) , and 
(12 — 3) X 4 = 36 if N G Had (12) . The reader will easily check that our theorem 
accurately predicts these values and appropriately assigns pn = 8 . 

When n = 8 , we have RS (8) = {0, 8} . In fact, a 3-normalized matrix G Had (8) 
will have exactly one row with sum 8 and seven with row sum 0. For, on the one 
hand, all eight rows cannot sum to zero; this would make them linearly dependent. 
On the other hand, all entries in a row with sum 8 must be +ls, so even two such 
rows create a dependency. We consequently have ex (N) = 8 , again confirming the 
theorem. □ 

Remark 4.3. According to Prop. 12.91 the excess of a 3-normalized matrix N € 
Had (n) always has certain congruence properties. The bound v* in Theorem 
14.21 however, may not have these properties, and in such cases, we cannot achieve 
ex (A^) = z/* . This occurs for the first time, for instance when n = 80 = (mod 16) . 
Then p„ = 12 and v* = 677^3- By Proposition 12.91 however, the excess of any 
actual 3-normalized N G Had (80) must be divisible by 16. Since the largest multi- 
ple of 16 below z^gQ is 672 , the latter provides a sharper upper bound for ex (A) 
when n = 80 . We know no 3-normalized A^ G Had (80) with excess above 624, 
but if there exists a matrix A^ G Had (80) with no = 1 (ie. 1 row with excess 0), 
ng = 68 and nig = 8 , it has excess 672, thus attaining the sharper bound. 

In general, when this situation arises (the next such occurs when n = 104), our 
bound can be sharpened similarly. The precise formulation, using \_x\ to denote 
the integer part of x is as follows: 



nripmin - r){pmm + 8 - r) = for all r G RS (n) . 



8 



- +4 if n 
2 



4 (mod 8) , 



ex (A) < < 




- +8 if n 
2 



8 (mod 16) 



V 




if n 



(mod 16) . 
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We should emphasize, however, that even this sharpened bound may not be attained 
for ah n . Indeed, it remains unknown whether Had (n) is non-empty for all n = 
(mod 4) . 

Comparison of bounds. As discussed at the beginning of this section, our con- 
struction (Theorem 13. 1|) produces a determinant larger than any the maximal ex- 
cess method can construct whenever there exists a 3-normalized Hadamard matrix 
A'^ G Had (n) for which (|24|) holds. Given the n^/^ bound on maximal excess, (|24|) 
certainly holds whenever 

(28) ex (A^) + n > n^^^ . 

We now show that Theorem 14.21 presents no obstruction to the existence of such 
a matrix A^ , except in the very lowest dimensions. Recall that v* , defined in 
that theorem, denotes our upper bound for the excess of a 3-normalized Hadamard 
matrix in Had (n) . 

Proposition 4.4. For all n > 12 we have 



(29) V(n-3)(n-4)(n- 12) < i^* < Wn - 3 , 

so that f* = n^/^ — 0{^/n) , and {i^* +n) — r?!'^ ~ n as n ^ oo in the sense that 

(30) lim ^ " ^ = 1 . 

In fact, for all n > 88 , we have 

(31) vl^n> n^/^ , 

Proof. It is clear from the definition of t'* (and that fact that pn> ^) that we can 
bound from below by the minimum of the function 

^ , , (n-3)x (n-4)(n-12) 

(32) Gn{x) = ^ ^ ' + ^ , X G (0, oo) . 

(We minimized essentially the same function in proving Theorem 14.21 but over a 
discrete domain.) 

When n > 12 , the function is concave up, with one critical point on (0, oo) . Com- 
puting the minimum explicitly on this interval, one easily obtains the lower bound 
in (|29|) when n > 12 . The bound is trivial when n = 12 . 



To get the upper bound z^* < n^Jn — 3 , we note that the definition of pn also 
implies 

+ 4 > p„ > ^ - 4 . 



\/n — 3 ~ ~ \/n — 3 
If we restrict x to this interval, the maximum of Gn{x) must occur at an endpoint, 
and an easy calculation shows that Gn{x) = n\/n — 3 at both endpoints. This 
value consequently bounds ly* from above, as claimed. 
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The asymptotic result u* = r?!'^ — 0{^/n) now follows. In fact, one easily uses the 
the bounds in (PU)) to compute 

3^1- K- n^''^ ^ . ,v*n- n3/2 19 
> nmsup -= — > iimmi — — > . 

2 n^+oo \n n^+oo Jn 2 



These facts make pO() obvious, and imply that holds for sufficiently large n . 

To show more precisely that H31|l holds whenever n > 88 , it will suffice to do so 
with f* replaced by the lower bound in (|29() : that is to show 

V(n - 3)(n - 4)(n - 12) > n^/^ _ ^ ^j. „ > gg . 

Squaring both sides and simplifying the result we find that this inequality is equiv- 
alent to 

- lOn^) + (48n - 72) > . 
When yjn > 10 , the first term in parentheses here is clearly non-negative, and the 
second is positive. We therefore see that (|31() holds for all n > 100 . By hand, one 
checks that it holds for n = 96 , 92 and 88 as well, but not for n = 84 . □ 

As mentioned earlier, Kounias & Farmakis bound the maximal excess of an order-n 
Hadamard by a number cj* < n^/^ when n is not a perfect square. For n < 88 , 
our bound v* + n dips below n'^/^ , but it still exceeds a* in most cases. Indeed, 
for all n = 4 (mod 8) , we have i'*+n>a^, with equality only for perfect square 
n . For n = (mod 8) , we have u* + n > a* for all n > 44 . For n < 100 we 
display the exact values of a* and i'*+n in Table 1 below. 



5. Examples 



Theorem 14 . 21 b ounds the excess of a 3- normalized Hadmard matrix from above, and 
by virtue of Theorem 13.11 this provides an upper bound for the determinants we 
can hope to achieve using the construction described at the beginning of 331 These 
facts would hold little interest if the construction did not actually produce examples 
that attain, or at least approach these upper bounds, thereby providing good lower 
bounds for the maximal determinant problem in the corresponding dimensions. 

Our construction does produce numerous examples of this type. With regard to low 
orders, it generates globally determinant-maximizing sign matrices of orders n = 5, 
13, and 21. By way of comparison, the maximal excess method does this for n = 5, 
9, and 17, but not for n = 13 or 21. For n = 45, 53, and 69, our method reproduces 
the largest known determinants, reported in jUSDSj . Most significantly, we have 
used it to set new records for n = 49, 65, 73, 77, 85, 93, and 97. 

Theorem 13. II makes it clear that to set records of this type using our construction, 
we need to find 3-normalized Hadamard matrices of large excess. In particular, to 
improve upon the maximal excess technique, we need to satisfy condition (|24l) . 
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The general literature of Hadamard matrices provides a variety of techniques for 
constructing Hadamard matrices, and explicit examples of relatively low order ( n < 
200 , say) abound on the internet. Unfortunately, we are unable to predict, based on 
construction technique or structural elements, for example, which of these matrices 
will have large excess after 3-normalization. 

We have therefore searched for large-excess 3-normalized Hadamard matrices using 
an approach with two main steps. Namely, 

(1) Generate initial Hadamard matrices from distinct equivalence classes in 
Had (n) , and then 

(2) Find a matrix in each class which attains the maximum excess subject to 
the 3-normalization constraint. 

A program we wrote to perform the latter task runs in polynomial time and has 
served us well for n < 100 or so. Finding initial Hadamards which can be 3- 
normalized with high excess, however, poses real difficulty in general. Starting in 
dimension 56 for the case n = , and dimension 92 for the case n = 4 (mod 8) , 
we have not yet succeeded in finding 3-normalized Hadamard matrices with excess 
achieving our upper bound . 

We next discuss our approach to the two tasks above. 



Generating initial Hadamards. For orders n = 4 (mod 8), online libraries of 
Williamson and so-called "good" matrices jSi] provide a rich source of initial Had- 
amard matrices. When n = (mod 8), however, one finds fewer examples online, 
and we have resorted to generating our own from smaller Hadamards. The tensor 
product of Hadamard matrices Hi £ Had (si) and H2 G Had (.52) always belongs to 
Had (S1S2) . More generally, one can construct Hadamard matrices of order S1S2 by 
a technique known as "weaving," given si matrices in Had (52) and S2 matrices in 
Had(si) (Crj . Craigen and Kharaghani have "woven" Hadamard matrices with very 
high excess |CKj . Unfortunately weaving does not, in general, produce 3-normalized 
Hadamard matrices of high excess. 

On the other hand, we have used another technique, the Multiplication Theo- 
rem of Agaian-Sarukhanyan, with quite good results. This Theorem constructs 
a Hadamard of order 8/01^2, given Hi G Had(4A;i) and H2 G Had (4^2)- The 
reader may consult "Kg or |SYj for further details. 



Excess-maximization. Our algorithm for maximizing excess on an equivalence 
class of 3-normalized Hadamard matrices relies on the fact that the excess depends 
only on a choice of the three rows used in the 3-normalization. We make this fact 
precise in Lemma 15.21 below, which requires the following definition: 
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Definition 5.1. Let denote the group of n x n signed permutation matrices, 
i.e. matrices obtained by permuting the rows of an n x n diagonal matrix of it Is. 
Let 5'^n_3 C denote the subgroup comprising matrices with (3, n — 3) block- 
diagonal form. 

Lemma 5.2. Suppose N is a 3-normalized Hadamard matrix of order n > 2 . 
Then any 3-normalized Hadamard matrix N' equivalent to N can be written as 

(33) N' = RNC, where R,CeS^, 

and ex ( A^') depends only on the coset to which R belongs in /5'^„_3 . 

Proof. Equation H33|) simply restates the definition of "equivalent" given in Defn. lT^ 
The main conclusion here is the one about excess, and for that, it suffices to show 
ex (RNC) = ex (A^) whenever both matrices are 3-normalized, and R £ S^^_^ . 

To do so, observe that rows 1,2, and 3 of form a 3 x n submatrix with columns 













+ , 


1 


+ 


, and 




+ 


_ + _ 









where +, — signify ±1 respectively. These vectors occur in successive blocks of 
length k := n/4, no two are mutually opposite, and signed row permutations pre- 
serve their distinctness. It follows that the signed permutation C in ()33() can 
permute columns within any one of the four /c-blocks, or move entire /c-blocks, but 
it must maintain the contiguity of each block. Further, one sees that C negates 
all columns or none, depending on whether R negates an odd or even number, 
respectively, of the first three rows. 

Now recall from Proposition 12.91 that for each i = 4, 5, . . . , n , row i oi N has the 
form (aj , bj , Cj , dj) , where aj, bj, Cj, dj G {±1}'^ all have equal excess. The facts 
above now force row i of A'^C to take the form 

i^i ) t)j , Cj , dj) , 

where we get {a^, b^, c'-, d^} by permuting {aj, bj, Cj, dj} , permuting the entries 
of each vector, and then, possibly, negating them all. Moreover, since R belongs 
to the subgroup S^j^_^ the last n — 3 rows of R ■ NC equal those NC , modulo 
permutations and sign changes of entire rows. It follows that for some permutation 
a of {4, 5, . . . , n} , we have 

(34) ex (a9 = ex (bQ = ex (c^ = ex (dQ = ±ex (a^(j)) , 

But the excesses of rows 4 through n of RNC must be non-negative fDefinition l2.8() . 
so the final sign in (|34|) must always be That is, rows 4 through n of N and 

RNC have exactly the same excesses, modulo permutation. These rows account for 
the entire excess of a 3-normalized Hadamard matrix, so we are done. □ 
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The subgroup S^^_^ C is clearly isomorphic to x 5^_3 , and the order of 
is clearly 2"n! for any n . We therefore have 

(35) # (s^/53±„_3) = (3) • 

Indeed, the (3) matrices we get by choosing three rows of the nxn identity matrix, 
then exchanging them (in any order) with rows 1, 2, and 3, each represent a different 
coset of S±/5j^_3 . 

Using these facts we easily find a 3-normalized matrix of maximal excess equivalent 
to any initial H G Had (n) . We take jj := (3) matrices {Ra}a=i C , each 
representing a different coset of / S^^_^ , and then form the matrices 

(36) RiH, R2H, R^H £ Had (n) 

According to Proposition I2.9[ we can 3-normalize each RiH using only a signed 
column permutation Cj , and some row-negations, which we can represent by a 
diagonal matrix Di of ±ls, multiplying Ri on the right. Such diagonal matrices 
clearly belong to 5'^„_3 , so we now have (3) 3-normalized Hadamard matrices Ni , 
i = 1,2,..., Q), with 

Ni = RiDiHCi . 

Since the RiDiS represent every coset of / Sf^_^ , Lemma ESI above guarantees 
that among these, the Ni with largest excess actually maximizes excess among all 
3-normalized Hadamard matrices equivalent to H . 

Remark 5.3. Already when n = 16 , one of the five equivalence classes in Had (16) 
fails to contain any 3-normalized matrix that achieves the bound in Theorem 14.21 

By tabulating the row sums one gets after 3-normalizing representatives from each 
of the (3) cosets described above, one associates to any H E Had (n) a statistic 
equivalent to the profile defined by Cooper, Milas, & Wallis in |CMWj . 

A table of large determinants. For each n = (mod 4) , 4 < n < 100 , Table 1 
below displays the results of our efforts to find 3-normalized Hadamard matrices of 
large excess, and the large determinants they generate via our construction. Addi- 
tionally, the table compares our results to those obtained using the maximal excess 
technique of Farmakis & Kounias, and indicates the largest known determinant of 
each order. The actual matrices which attain these determinants are posted on our 
website |OSj . 

We now give a detailed explanation of Table 1. 
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Upper bounds 


Best known 








n + 1 


n + u* 


< 










fJ-n+l/Pn+l 


5 


8 


8 


8 


8 


3* 


3 


1 


9 


16 


20 


16 


20 




8.25 


0.85 


13 


48 


36 


48 


36 


15* 


15 


1 


17 


48 


64 


48 


64 


20* 


22.98 


0.87 


21 


96 


80 


96 


80 


29* 


32.02 


0.91 


25 


96 


112 


96 


112 


42* 


42 


1 


29 


152 


140 


152 


140 


4577 


52.85 


0.87 


66 


16U 


172 


loO 


172 


55 /s 


64.50 


0.85 


37 
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216 


216 


72 
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76.90 


0.94 
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Ub 


1 1 Q 1 O 

iio.iy 
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obo 


od4 


obo 
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lUO 


1 QQ 01 

ioo.zi 
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( 


>l >l O 

44o 
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oo4 


A nn 
4UU 


11/1 
ii4 


1 /I Q QO 
148. oZ 


n 77 

U. M 


Di 


456 


452 


456 


440 


iOO 


1 

iUO 


1 


65 


560 


512 


528 


512 


148 


181.73 


0.81 


69 


552 


548 


552 


544 


155 


198.98 


0.78 


73 


656 


600 


624 


580 


174 


216.75 


0.80 


77 


656 


652 


656 


628 


183 


235.02 


0.78 


81 


752 


704 


704 


704 


196 


253.77 


0.78 


85 


768 


756 


768 


756 


213 


273 


0.78 


89 


864 


812 


768 


792 


220 


292.69 


0.75 


93 


888 


872 


864 


828 


239 


312.83 


0.76 


97 


976 


932 


928 


920 


256 


333.42 


0.77 


101 


1016 


1000 


984 


1000 


275 


354.44 


0.78 



Table 1. Largest known excesses of order n = (mod 4) and de- 
terminants of order n + 1 . //n+i encodes the largest known determi- 
nant of order n + 1 . Compare Hn+i with the corresponding factor 
for our determinant ( ^"+7^+" ), and that of [FKI ( ^ ). 



(T* and an respectively denote the sharpest known upper bound for the 
maximal excess, and the largest known excess, of a Hadamard matrix of 
order n . 

f* and f„ respectively denote the sharpest known upper bound for the 
excess of, and the largest known excess for, a 3-normalized Hadamard matrix 
of order n . 

fin+i and (3n+i abbreviate the largest known determinant detM^+i, and 
the Barba bound B{n + 1) (see (0)), by dividing out large integer factors 
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as follows: 

_ det Mn+i _ B{n + l) 

Here k = n/A, and M„+i is a sign matrix of order n + 1 having largest 
known determinant. 



Explanatory notes. 

(1) The "Upper hounds'' columns. The tabulated values for c* all come from 
|l^'Kl Table 2]. Except for n = 80 (see Remark |4.3() . we computed the values 
for J/* using Theorem l4.2[ and we list n + i^* instead of just z^* , to facilitate 
the key comparison indicated by In each row, boldface indicates the 
larger of re + t'^ and u* . 

(2) The "Best known" columns. A bold entry in either column indicates that 
the largest known determinant of that order is produced by applying the 
corresponding construction to a matrix with the associated excess Vn or 
Gn . When neither entry is bold, neither method constructs the current 
record. 

All values for o"„ come from |FK1 Table 2] except for a^2 = 580 and a-j^ = 
628 , which were reported in |OSDSj . 

(3) The "Hn+i " column. Current records which can be constructed by our 
method are underlined. A bold entry represents a determinant that achieves 
a new record, reported for the first time here. Values of fin+i corresponding 
to determinants known to attain md(n + 1) are marked with an asterisk. 

(4) For n < 20 , either the 3-normalized or maximal excess method (or both) 
constructs an order n + 1 sign matrix that achieves the global determinant 
maximum. Equivalent matrices were all constructed by earlier investigators, 
however: See [MollWil for n + 1 = 5 . jET) for n+1 = 9 . for n + 1 = 13 . 
jMEl for n + 1 = 17 , and \CKM\ for rt + 1 = 21 . 

(5) For sizes n + 1 = 25, 29, 33, 37, 41 and 61 , neither our method nor that of 
Farmakis & Kounias constructs the largest known determinant. In dimen- 
sions 25,41, and 61, matrices coming from block designs (SBIBDs) attain 
100% of Barba's bound B{n + 1) ; see jBH Eil for n + 1 = 25 , BHHj for 
n + 1 = 41 , and |Brj for n + 1 = 61 . The largest known determinant for 
n+1 = 33 was reported in |OSDSj . For n + 1 = 29 and 37, we sketch a 
construction for the largest known determinants below. 
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6. Improved results for n + 1 = 29, 37 

The "Upper bounds" and "Best known" columns of Table 1 show that for n+1 = 29, 
our construction achieves a determinant larger than any the maximal excess method 
can produce. The value /i29 = 45^/7 there, however, corresponds to the still larger 
determinant reported in |()SDSj using a gradient ascent algorithm. We have since 
discovered that the latter matrix can be formed by applying a symmetric version 
of our construction. To do this, we start with a matrix H E Had (28) for which 
both H and satisfy a generalized 3-normalization condition, and for which 
both have the optimal set of row sums prescribed by Theorem 14.21 We then apply 
appropriately adapted versions of steps 1-4 of our construction to H , transpose 
the result, and apply the same four steps again. This produces a matrix equivalent 
to the one reported in ^OSDS^. 

Hoping to profit by applying the same idea to other orders, we have so far achieved 
one notable success: A matrix A £ {±1}^'^^^'' with det^ = 72 x 9^'^ x 2^^. This 
determinant, reported here for the first time, attains 94% of the Barba bound, and 
we conjecture it to be the global maximum for its order. By way of comparison. 
Table 1 shows that the largest determinant that either the maximal excess method 
or our own can hope produce is 63 x 9^^ x 2^^. This value has been achieved by 
both methods, but it represents just 82% of Barba's bound. 

Guided by the n+1 = 29 case, we start by finding a matrix H £ Had (36) which can 
be "doubly 3-normalized" with optimal row and column sums. In this dimension, 
the three 3-normalizing rows and columns overlap in a somewhat more complicated 
way than in dimension 28. This forces further modifications of the double 4-step 
construction described above, but the procedure is very similar. 

We are currently investigating the relationship between these two examples and 
possible generalizations. We will publish further details of this work in our sequel to 
the present work. Meanwhile, we have posted both the n + 1 = 29 and n + 1 = 37 
examples on our website |OSj . 
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